Structure of propagators for quantum nondemolition systems 
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In the scheme of a quantum nondemolition (QND) measurement, an observable is measured 
without perturbing its evolution. In the context of studies of decoherence in quantum computing, 
we examine the 'open' quantum system of a two-level atom, or equivalently, a spin-1/2 system, 
in interaction with quantum reservoirs of either oscillators or spins, under the QND condition of 
the Hamiltonian of the system commuting with the system-reservoir interaction. The propagators 
• for these QND Hamiltonians are shown to be connected to the squeezing and rotation operators 

for the two baths, respectively. Squeezing and rotation being both phase space area-preserving 
canonical transformations, this brings out an interesting analogy between the energy-preserving 
QND Hamiltonians and the homogeneous linear canonical transformations. 
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In the scheme of a quantum nondemolition (QND) measurement, an observable is measured without perturbing its 
free motion. Such a scheme was originally introduced in the context of the detection of gravitational waves It was 
to counter the quantum mechanical unpredictability that in general would disturb the system being measured. The 
dynamical evolution of a system immediately following a measurement limits the class of observables that may be 
measured repeatedly with arbitrary precision, with the influence of the measurement apparatus on the system being 
confined strictly to the conjugate observables. Observables having this feature are called QND or back-action evasion 
observables @,M 0]- ^ n addition to its relevance in ultrasensitive measurements, a QND scheme provides a way to 
' prepare quantum mechanical states which may otherwise be difficult to create, such as Fock states with a specific 
y—i , number of particles. One of the original proposals for a quantum optical QND scheme was that involving the Kerr 
medium [f|, which changes its refractive index as a function of the number of photons in the 'signal' pump laser. 
The advent of experimental methods for producing Bose-Einstein condensation (BEC) enables us to make progress 
■ in the matter-wave analogue of the optical QND experiments. In the context of research into BEC, QND schemes 
fS |' with atoms are particularly valuable, for instance, in engineering entangled states or Schrodinger's cat states. A state 
! 1 1 , preparation with BEC has recently been performed in the form of squeezed state creation in an optical lattice [||. 

In a different context, it has been shown that the accuracy of atomic interferometry can be improved by using QND 
& , measurements of the atomic populations at the inputs to the interferometer Q. 

No system of interest, except the entire universe, can be thought of as an isolated system - all subsets of the 
universe are in fact 'open' systems, each surrounded by a larger system constituting its environment. The theory of 
open quantum systems provides a natural route for reconciliation of dissipation and decoherence with the process of 
quantization. In this picture, friction or damping comes about by the transfer of energy from the 'small' system (the 
system of interest) to the 'large' environment. The energy, once transferred, disappears into the environment and is 
not given back within any time of physical relevance. Ford, Kac and Mazur Q suggested the first microscopic model 
describing dissipative effects in which the system was assumed to be coupled to a reservoir of an infinite number of 
harmonic oscillators. Interest in quantum dissipation, using the system-environment approach, was intensified by the 
works of Caldeira and Leggett and Zurek [To| among others. The path- integral approach, developed by Feynman 
and Vernon 11 1, was used by Caldeira and Leggett [9(], and the reduced dynamics of the system of interest was 
followed taking into account the influence of its environment, quantified by the influence functional. In the model of 
the fluctuating or "Brownian" motion of a quantum particle studied by Caldeira and Leggett the coordinate of 
the particle was coupled linearly to the harmonic oscillator reservoir, and it was also assumed that the system and 
the environment were initially factorized. The treatment of the quantum Brownian motion has since been generalized 
to the physically reasonable initial condition of a mixed state of the system and its environment by Hakim and 
Ambegaokar [12| . Smith and Caldeira [13j| . Grabert, Schramm and Ingold [3], and by us for the case of a system in a 
Stern-Gerlach potential [l| , and also for the quantum Brownian motion with nonlinear system-environment couplings 
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An open system Hamiltonian is of the QND type if the Hamiltonian Hs of the system commutes with the Hamilto- 
nian Hsr describing the system-reservoir interaction, i.e., Hsr is a constant of motion generated by H$ . Interestingly, 
such a system may still undergo decoherence or dephasing without any dissipation of energy (l7l fl8l]. 

In this Letter, we study such QND 'open system' Hamiltonians of particular interest in the context of decoherence 
in quantum computing, and obtain the propagators of the composite systems explicitly using path integral methods. 
The aim is to shed some light on the problem of QND measurement schemes. Can one draw upon any familiar 
symmetries to connect with the time-evolution operation of these QND systems? 

We take our system to be a two-level atom, or equivalently, a spin-1/2 system. We consider two types of environment, 
describable as baths of either oscillators or spins. One cannot in general map a spin-bath to an oscillator-bath (or 
vice versa) ; they constitute distinct universality classes of quantum environment [1 91 ] . The first case of oscillator-bath 



models (originated by Feynman and Vernon 11[) describes delocalized environmental modes. For the spin-bath, on 
the other hand, the finite Hilbert space of each spin makes it appropriate for describing the low-energy dynamics 
of a set of localized environmental modes. A difficulty associated with handling path integrals for spins comes from 
the discrete matrix nature of the spin-Hamiltonians. This difficulty is overcome by bosonizing the Hamiltonian by 
representing the spin angular momentum operators in terms of boson operators following Schwinger's theory of angular 
momentum [20j. 

We then use the Bargmann representation [2l[ for all the boson operators. The Schrodinger representation of 
quantum states diagonalizes the position operator, expressing pure states as wave functions, whereas the Bargmann 
representation diagonalizes the creation operator b\ and expresses each state vector \tp) in the Hilbert state H as an 
entire analytic function f(a) of a complex variable a. The association — > /(<*) can be written conveniently in 
terms of the normalized coherent states |a) which are the right eigenstates of the annihilation operator b: 

b\a) — a\a), 

(a'\a) = exp (-^\a'\ 2 - ^\a\ 2 + a'*a 



giving 



We obtain the explicit propagators for these many-body systems from those of the expanded bosonized forms by 
appropriate projection. 

The propagators for the QND Hamiltonians with an oscillator bath and a spin bath are shown to be connected to the 
squeezing and rotation operators, respectively, which are both phase space area-preserving canonical transformations. 
This brings out an interesting analogy between the energy-preserving QND Hamiltonians and the homogeneous linear 
canonical transformations, which would need further systematic probing. 

We first take the case where the system is a two-level atom interacting with a bosonic bath of harmonic oscillators 
with a QND type of coupling. Such a model has been studied [12, 0, HiJ in the context of the influence of decoherence 
in quantum computation. The total system evolves under the Hamiltonian, 

H = H s + H R + H S r 

= ^r°"* + £ M{b k + (^) Y,9k{h + bl)a z . (1) 

fc=l ^ ' k=l 

Here Hs,Hr and Hsr stand for the Hamiltonians of the system, reservoir, and system-reservoir interaction, respec- 
tively. We have made use of the equivalence of a two-level atom and a spin-1/2 system, a x , a z denote the standard Pauli 
spin matrices and are related to the spin- flipping (or atomic raising and lowering) operators S+ and S-: a x = S+ + S-, 
a z — 25+ S- — 1. In {J) b k ,bk denote the Bose creation and annihilation operators for the M oscillators of frequency 
u>k representing the reservoir, stands for the coupling constant (assumed real) for the interaction of the field with 
the spin. Since [Hs,Hsr] = 0, the Hamiltonian |T]) is of QND type. 

The explicit propagator exp(— ^p) for the Hamiltonian JT]) is obtained by using functional integration and bosoniza- 
tion [25l.l26l|. the details of which can be found elsewhere [27] . In order to express the spin angular momentum operators 
in terms of boson operators, we employ Schwinger's theory of angular momentum [20( by which any angular momen- 
tum can be represented in terms of a pair of boson operators with the usual commutation rules. The spin operators 
a z and a x can be written in terms of the boson operators ag, and a 7 , a 7 as 

a z = at,a 7 - a) ap, 
t t 
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In the Bargmann representation [2l[ the actions of b and are 

&t/(«) = a *f(a), 



- ( 2 ) 



where |a) is the normalized coherent state. The spin operator becomes 

Here the variable (3* is associated with the spin-down state and the variable 7* with the spin-up state. 

The propagator for the bosonized form of ([T]) is given as a path integral over coherent state variables, and the 
required propagator is then obtained by appropriate projection, the amplitudes of which are given in a matrix form 
as [13 



Ul (a*,/r, 7 V;a',/?',y,0) = exp V^^e^*' 




(4) 

In the above equation a is a vector with components {otk\ and 

\2J k=i uj k \ZJ k=i uj k 

M 

B = J2<Pk(a* k + a' k )+i^t, (6) 

k=l 

<fe = ^(l-e-^). (7) 

Here we associate the values a* with time t and a 1 with time t — 0. The simple form of the last term on the right-hand 
side of Q reveals the QND nature of the system-reservoir coupling. Since we are considering the unitary dynamics of 
the complete Hamiltonian fl} there is no decoherence, and the propagator ([!]) does not have any off-diagonal terms. 
In a treatment of the system alone, i.e., an open system analysis of Eq. (fTJ) after the tracing over the reservoir degrees 
of freedom, it has been shown [23| that the population, i.e., the diagonal elements of the reduced density matrix of 
the system remain constant in time while the off-diagonal elements that are a signature of the quantum coherences 
decay due to decoherence, as expected. 

Next we consider a variant of the Hamiltonian (fl]) wherein we include an external mode in resonance with the 
atomic transition, such as 

H = —^-a z + hila [ a 2~ az 

m /hu\ M 

+ Y J Mlb k +[ — \Y J 9k{hk + b k )a z . (8) 

k=l ^ ' k=l 

Here 

Cl = 2e.d*, (9) 

where d is the dipole transition matrix element and e comes from the field strength of the external driving mode E^lt) 
such that 

E L {t) = ee- Mt + e*e luJt . (10) 
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The use of <x z instead of the usual a x in the third term on the right-hand side of ([8]) allows us to have a QND 
Hamiltonian. Proceeding as before and introducing the symbol v* for the external mode we obtain the amplitudes 
of the propagator for fH} as 



U-2 



expj.V'e-^} e A ( e l 2 e \ ) , (11) 



where A is as in Eq. ([5]), 



B 2 = J2<f>k(<x*k+cx'k)+i (—z -)t, (12) 
fe=i ^ ' 

4>k is as in Eq. 0. 

Now we consider the case where the reservoir is composed of spins or two-level systems, as has been dealt with by 
Shao and collaborators in the context of QND systems [18| and also quantum computation [28| , and for a nanomagnet 
coupled to nuclear and paramagnetic spins [191 ]. The total Hamiltonian is taken as 



H — Hs + He + Hsu 

huj M hui M 

= —S z + s ^Tiuj k (T zk + — s ^c k (j xk S z . (13) 

k=l k=l 



Here we use S z for the system and a zk , a xk for the bath. Since [Hg, Hsr] = 0, we have a QND Hamiltonian. In the 
Bargmann representation, we associate the variable 3* with the spin-down state and the variable 7* with the spin-up 
state for the bath variables, and we have 

— > * 9 3* 9 
dj* dd* ' 

°-^+W + rw- (14) 

Similarly, the bosonization of the system variable gives 

s-^e£-rw> as) 

where the variable 8* is associated with the spin-down state and the variable £* with the spin-up state. The amplitudes 
of the propagator for (|13p . in the Hilbert space of Hr, are obtained as 



m 00 r T ? T f 



k=l n=0 







Vp <fS.t/ cos9 fe (") zsin6 fc (") , 

' ' I (-l^sinG^ (-l)"cose fc (™) ' ' 



Here f3* , 7* are vectors with components {8 k } and {7*:}, respectively, and 



e *(») = _ Sz c k A ni (17) 



^ = E(-l) J+1 2r, + (-irt. (18) 
3=1 

In Eq. (|16[) . expanding the terms containing S z in the system space, we can see that out of the 16 amplitudes of the 
propagator for each mode of the reservoir, only the energy-conserving diagonal terms survive in Eq. (|16|) due to the 
QND nature of the system-reservoir coupling. 
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We look closely at the forms of the propagators (|4]) and (p~6|) of the QND type Hamiltonians |T]) and (| 13|) , respectively. 
In the first case with an oscillator bath, Eq. (JD) involves the matrix 

e B 
e~ B 

where B is given by Eq. This can be used to generate the following transformation in phase space: 

X\ ( e B \ ( x 



PI \ e- B I \p 



(19) 



It can be easily seen from Eq. (fT9|) that the Jacobian of the transformation is unity and it is a phase space area- 
preserving transformation. The first matrix on the right-hand side of (|19p has the form of a 'squeezing' operation 
[29j, which is an area-preserving (in phase space) canonical transformation coming out as an artifact of homogeneous 
linear canonical transformations [3Cj • 

In the second case of a spin bath, Eq. ([To]) involves the matrix 

/ cose fe (") isin9 fe (™) \ , s 

U - ^(-l) n isine fc (") (-l) n oose fc C»»)J ' 

where 6 fc( ") is given by Eq. (fT"7|l . For particular n and k, we write O fe ( n ) as 8. For n even, the above matrix (|20|) 
becomes 



cos 6 i sin 6 
i sin 9 cos 6 



Using the Campbell-Baker-Hausdorff identity [3l| this matrix can be shown to transform the spin vector a = 
(a x ,a v ,a z ) as 

/I \ /CT a \ 

= cos26 -sin26 \ [ cr y , (22) 

\0 sin 26 cos 26 / \a z J 

i.e., the abstract spin vector is 'rotated' about the x-axis by an angle 26. For n odd, (|2U)l becomes (again writing 
@ fc («) f or particular n and k as 6) 




cos 6 i sin 6 \ _ / cos 6 i sin 6 
—i sin 6 — cos 6 / 2 \ i sin 6 cos 6 



a z e ie ^. (23) 



Thus the n-odd matrix is related to the n-even matrix by the spin-flipping energy. The above matrix transforms the 
spin vector a as 

rr z e ie ^ ( a y ) e~ i& ^a z = e" (o cos 26 sin 26 \ ( a y \ (24) 

VO sin 26 -cos 26/ \a z J 

It can be easily seen from the right-hand side of the Eq. that the determinant of the transformation of the 

spin vectors brought about by the n-odd matrix (|23[) has the value unity. It is well known that the determinant of a 
rotation matrix is unity [32j • Thus we see that the above transformation has the form of a rotation. Specifically, it 
can be seen that 




cos 26 sin 26 \ _ ( cos 26 sin 26 
sin 26 -cos 26 J ~ I - sin 26 cos 26 



(25) 



and 

/ cos 26 sin26\ T /cos 26 -sin 26 \ , , 

-sin 26 cos 26 J _ I sin 26 cos 26 J' [ ' 



Here T stands for the transpose operation. From the above it is seen that the matrix (I2TH) has the form of the operation 
of 'rotation', which is also a phase space area-preserving canonical transformation [29j and comes out as an artifact 
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of homogeneous linear canonical transformations [3(j. Any element of the group of ho mog eneous linear canonical 
transformations can be written as a product of a unitary and a positive transformation |33l |34| , which in turn can 
be shown to have unitary representations (in the Fock space) of rotation and squeezing operations, respectively [30j ] . 
It is interesting that the propagators for the Hamiltonians given by Eqs. |T]) and (fT3")l . one involving a two- level 
system coupled to a bath of harmonic oscillators and the other with a bath of two-level systems, are analogous to the 
squeezing and rotation operations, respectively. 

In conclusion, we have investigated the forms of the propagators of some QND Hamiltonians commonly used in 
the literature, for example, for the study of decoherence in quantum computers. We have calculated the propagators 
for the Hamiltonian of a two-level system interacting with a bosonic bath of harmonic oscillators (and its variant 
wherein we have included an external mode in resonance with the atomic transition), as also a spin bath of two- level 
systems. In each case the system-bath interaction is of QND type, i.e., the Hamiltonian of the system commutes with 
the Hamiltonian describing the system-bath interaction. We have found an interesting analogue of the propagators 
of these many-body Hamiltonians to squeezing and to rotation, respectively. Every homogeneous linear canonical 
transformation can be factored into the rotation and squeezing operations and these cannot in general be mapped 
from one to the other - just as one cannot in general map a spin bath to an oscillator bath (or vice versa) - but together 
they span the class of homogeneous linear canonical transformations and are 'universal'. Squeezing and rotation, being 
artifacts of homogeneous linear canonical transformations, are both phase-space area-preserving transformations, and 
thus this implies a curious analogy between the energy-preserving QND Hamiltonians and the homogeneous linear 
canonical transformations. This insight into the structure of the QND systems would hopefully lead to future studies 
into this domain. 
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